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The Comment by [T[ on our paper Q is essentially 
based on two claims: (i) The integrals used in [2J to de- 
termine the QNMs do not exist for negative values of the 
argument, thereby invalidating the extraction of poles 
using Born approximation, (ii) The poles of the integral 
used in Q to determine the quasinormal modes (QNMs) 
in the Schwarzschild case come not only from Gamma 
function but from Whittakcr functions as well [eq (7) of 

The second claim is easy to dispose of as erroneous, 
which we will do first. It is straightforward to see from 
the formula [9.211.4] of that 
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where 'J is the confluent hypergeometric (Tricomi) func- 
tion [also written as U(a, b, z)]. This function \fr is regular 
for all finite a and b. Hence, the only poles arise from 
the Gamma function and no other functions are needed 
for determining the pole structure. 

Let us now take up the first point (i) which is essen- 
tially related to the existence of integrals of the form [e.g., 
eq (12) of Q] 



dx x ^/*e luJx/K 



(2) 



where k = (4M)" 1 . If we introduce the parameters v — 
1 + ilj I ' k and // = —iuj/k, the integral can be evaluated 
as 



dx x v - x ^ x 



Re(n) > 0, Re(v) > 
(3) 

The condition Re(/i) > translates to ujj > 0, which 
does not affect our results as we are interested only in 
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(large) positive values of loj. If we let z = [xx, then we 
essentially have to evaluate the Gamma function integral 



H V I = I 1 



dz z v 



(4) 



The authors claim that this evaluation of the integral 
as Gamma function is valid only for Rev > (which 
translates to translates to u>i < k) while we are interested 
in ref. Q for large wj. This objection, too, is easy to take 
care of. 

The point to note is that, T(y) can be defined by ana- 
lytic continuation for negative values of v. Even though 
for Re(i^) < 0, the integrand in Eqf4] behaves as ~ z v ~ x 
as z — ¥ the analytic continuation allows one to define 
the Gamma function everywhere in complex plane. Once 
this is done, one can easily extract the poles. A simple, 
text book way to extract this result is to use the identity: 



i» r(-iz) = - 



(5) 



which provides a way to analytically continue F(i/) to 
negative values of Re(i^). This also shows that the func- 
tion has simple poles for negative integral values of Re(^) 
arising from the sin(7w) factor, which is also a well-known 
result. These are precisely the poles that we are interested 
in which gives us the desired QNM structure. Of course, 
by the very definition of a 'pole', the integral diverges 
at the pole; we stress that the whole exercise is to de- 
termine precisely where this occurs! The integral exists 
in a open neighborhood of the first order pole which is 
what we used in our analysis. The analytic continuation 
is based on the standard assumption in scattering theory 
that the scattering amplitude (and hence the integral in 
Eq. [5]) is analytic everywhere in the u> plane, except for 
a finite number of poles. The key idea developed in Q 
and Q and used in the paper under discussion Q was 
to use this assumption, identify the poles and relate it to 
the QNM. 

If one does not want to use the identity in Eqj5] but 
want to work directly with integral in EqU and extract 
the information about the poles, that is also possible. 
We only have to treat the integral in EqH] as a limit 
of a sequence of integrals with a suitable regularization 
parameter and study the poles. This can be done in 
many ways and we outline one procedure: Consider the 
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I' = / dz z"- 1 e- z e- a / z = 2a v l 2 K_ v ( c lJa\) (6) 
Jo 

which is well defined even for Re(^) < because of the 
regulator factor e^~ a ^ z \ (Here K v is the modified Bessel 
function and the relation can be obtained from Sec 
8.40-8.43.) We treat the integral in Eq|4] (especially for 
Re(^) < 0) as the limit of V when a is a positive infinites- 
imal quantity. This gives (again using Sec 8.40-8.43 of [3j 
and interpreting a — > + as a positive infinitesimal value 
for the regulator): 

poo 

I' = lira / dz z v ~ x q- z e~ a/z 

= lim 2a v/2 KJ2y/a) 

= I» + lim a v [T(-v) + 0(a)] (7) 

a->0+ 

Note that the last equality in the above equation is valid 
only when v is not a integer. In fact, if we take the limit 
v — > —n (where n is positive integer), the integral di- 
verges as it should, because our previous analysis using 
Eqj5]has already told us that the integral has simple poles 
at v = —n. We can determine the nature of the singu- 
larity arising from these poles trivially. When v — > — n, 
we obtain 

lim (y + n)I' = lim {v + n)T(i') 

v—t—n v^r — n 



+ lim lim a? [{v + n)T{-v) + 0{a)] 

a^0+ v— n 

= + lim a-' 1 ... ^ (8) 



which shows that the singularity of I 1 at v — — n is a 
simple pole for finite regulator. (The procedure is very 
similar to the ie prescription used in field theoretic calcu- 
lations.) Same arguments are valid for the integral (18) 
of We stress that the integrals are not expected to 
exist for v = —n, which are precisely the poles we want 
to determine! What we need is a sensible definition of 
the integral in the open neighborhood of the poles - 
which can be provided in many ways, of which we have 
described two. 

Finally, we would like to point out that paper in ques- 
tion which is being commented upon Q is a follow-up of 
two earlier papers Q and [f| developing the same tech- 
nique and containing the same integrals. It is somewhat 
surprising that the authors of [lj decided to comment 
a third, follow-up paper rather than the first two! We 
did not discuss the details of the regularization, analytic 
continuation etc. in our work [2j as the basic ideas were 
already implicit in the previous published work. 
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